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Necessary and sufficient conditions for bipartite entanglement are derived, which apply to arbi- 
trary Hilbert spaces. Motivated by the concept of witnesses, optimized entanglement inequalities 
are formulated solely in terms of arbitrary Hermitian operators, which makes them useful for ap- 
plications in experiments. The needed optimization procedure is based on a separability eigenvalue 
problem, whose analytical solutions are derived for a special class of projection operators. For gen- 
eral Hermitian operators, a numerical implementation of entanglement tests is proposed. It is also 
shown how to identify bound entangled states with positive partial transposition. 
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I. INTRODUCTION 

Entanglement is considered to be the key resource of 
the rapidly developing fields of Quantum Information 
Processing and Quantum Computation, for an introduc- 
tion see e.g. [ij- Among the early proposals concerning 
the applications of entangled states are those for quan- 
tum key distribution ,2], quantum dense coding Q , and 
quantum teleportation [J| . Despite the resulting great in- 
terest in entangled quantum states, their complete char- 
acterization is still an unsolved problem. 

It is known that entanglement can be fully identi- 
fied by applying all positive but not completely positive 
(PNCP) maps to a given state Q. The problem of this 
approach, however, consists in the fact that the general 
form of the PNCP maps is essentially unknown. The 
presently best studied PNCP map is the partial trans- 
position (PT) 0. It is known that PT gives a com- 
plete characterization of entanglement in Hilbert spaces 
of dimension 2 ® 2 and 2 ® 3 |5:]. Bipartite entangle- 
ment can also be completely characterized via PT in 
infinite-dimensional Hilbert spaces, as long as only Gaus- 
sian states are considered, whose moments up to sec- 
ond order fully describe their properties 0, [l] ■ By using 
higher-order moments, a complete characterization has 
been given for those entangled states which exhibit neg- 
ativities after application of the PT map [9|. This ap- 
proach gives no insight into bound entangled states re- 
maining positive after PT. To overcome this limitation, 
to the matrices of moments other kinds of PCNP maps 
have been applied llOl , in cluding Kossakowsky, Choi and 
Breuer maps [ll|, 13 ■ The identification of bound 
entanglement in this way, however, turned out to be a 
cumbersome problem. 

An equivalent approach of identifying entanglement is 
based on special types of Hermitian operators, called en- 
tanglement witnesses. The witnesses were introduced as 
a class of linear operators, whose mean values are non- 
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negative for separable states but can become negative 
for entangled states 0]. Presently only some classes of 
entanglement witnesses are available. Once a witness is 
known, an optimization can be performed 11411 . Also non- 
linear witnesses have been studied [H, [H, Il7( . which may 
increase the number of entangled states to be identified 
by a single witness, in comparison with a given linear 
witness. However, if one is able to construct the general 
form of the linear witnesses, the problem of identifying 
entanglement is completely solved. 

In the present contribution we show that any entan- 
glement witness can be expressed in terms of completely 
positive Hermitian operators, whose general form is well 
known. On this basis we derive entanglement inequali- 
ties, which are formulated solely in terms of general Her- 
mitian operators. We also provide an approach for op- 
timizing such inequalities, by introducing a separability 
eigenvalue problem. Our method is a powerful tool for 
analyzing experimental data, to verify any kind of entan- 
glement. One may also identify general bound entangled 
states whose density operator has a positive partial trans- 
position (PPT). 

The paper is structured as follows. In Sec. |TT] we de- 
rive the most general form of entanglement conditions in 
terms of Hermitian operators. This leads us to an opti- 
mization problem - the separability eigenvalue problem 
- which is studied in Sec. IIIII In Sec. IIVI this problem 
is solved for a class of projection operators and a gen- 
eral numerical implementation of entanglement tests for 
arbitrary quantum states is given. A method for the 
identification of bound entangled states is considered in 
Sec. |Vl Finally, in Sec. IVII we give a brief summary and 
some conclusions. 



II. ENTANGLEMENT CONDITIONS 

Let us consider two systems A and B, represented by 
arbitrary Hilbert spaces TLa and TLb with orthonormal 
bases being {\ei)}^^j and respectively, with / 

and J being arbitrary sets. Note that the Hilbert spaces 
are not necessarily finite or separable. Even spaces with 
an uncountable bases are under study. 
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An entanglement witness is a bounded Hermitian op- 
erator W, which has positive expectation values for sep- 
arable states and it has negative eigenvalues [1]. For our 
purposes a generalization of the class of entanglement 
witnesses is useful. One can think of bounded Hermitian 
operators, which have positive expectation values for sep- 
arable states: 

tr(o-W^) > (V(T separable), (1) 
||W^|| = sup{|w| : w eigenvalue of W} < oo. (2) 

All operators fulfilling the conditions ^ and ([2|) shall 
define the set VPab, operators in this set are called 
partial positive operators. All operators C fulfilling the 
conditions (^IC*!*) > (V|*) e Ha Hb) together 
with Eq. ([2]), with C in place of W, shall denote the 
set V of positive semi-definite operators. So all entan- 
glement witnesses are elements of the difference of these 
sets: VVab\V. 

It was shown by the Horodeckis [5l|, that for any en- 
tangled state Q there exists an entanglement witness 
W e VVab\P, so that the expectation value becomes 
negative: tr(^W^) < 0. For this inseparability theorem 
only linear entanglement witnesses were used, which are 
sufficient to identify all entangled states. For this reason 
we restrict our considerations to linear witnesses, which 
are elements of the set VPab- 

Let us consider the important example of witnesses 
based on the FT map. Recently it has been shown 9], 
that for any state g with a negative FT (NFT) there 
exists an operator /, such that 

{{PfD = tmPfD < 0. (3) 

These operators have been studied in detail as func- 
tions of the annihilation and creation operators of two 
harmonic oscillators: f{a,a\b,w). All the resulting 
iP f)^'^ are examples of elements of Wab, in particular 
they include all entanglement witnesses for NFT states. 

Now we will turn to the construction of entanglement 
witnesses in their most general form. As outlined above, 
the problem of finding all entanglement witnesses via 
FNCF maps is very difficult. Here we will introduce a 
different but equivalent approach to entanglement wit- 
nesses, which requires the class of V operators only. A 
Hermitian operator C is positive, if and only if it can be 
written as 

C = PI (4) 

In the first step we will now generate any entanglement 
witness out of a difference of positive operators. 

Lemma 1 For any entanglement witness W exists a real 
number A > and a positive Hermitian operator C so 
that W can be written as W — XI — C . 

Proof. The bounded operator in spectral decomposition 
is W — w dP{w), with P being a projection- 

valued measure and a{W) the bounded set of eigen- 
values. Let the supremum of all eigenvalues be 



A = supcr(W^). For all separable quantum states, 
W must be a positive map, which implies A > 0. 

W = X I 1 dP{w) - / (A - w) dP{w) 

=c, c>0 

= Ai - /" c dP{c) =\1-C, 

which is the demanded form. 

To formulate a new entanglement theorem for positive 
Hermitian operators, we need the definition of optimal 
entanglement witnesses as given by Lewenstein et al. [l^ . 
An entanglement witness Wi is finer than W2 , if and only 
if the entanglement of any state detected by W2 is also 
detected by Wi (beside other entangled states, which are 
not detected by W2). An entanglement witness W^opt is 
optimal, if and only if no witness is finer than Wopt- 

Therefore a state is separable, if and only if for all 
optimal entanglement witnesses the expectation value is 
positive. To find these witnesses, we need the function 
/ab(^), which maps a general Hermitian operator A to 
its maximal expectation value for a separable state: 

/ab(^) =sup{tr((TA) : a separable} 

==sup{(a,6|i|a,6) : {a\a) = {b\b) = 1}. (5) 

It is obvious, that a Hermitian operator VF = AI — (7 is 
a general element of PPab, if and only if A > fAB{C)- 
And it is optimal, if and only if A = fAB{C). 

Theorem 1 A state g is entangled, if and only if there 
exists C e P: fAB{C) < tr{gC). 

Proof. Let Wopt be an optimal witness, which detects the 
entanglement of g: 

trigWopt) < 
tY{gifAB{C)i-C))<0 

fAB{CMg)<tr{gC) 
fAB{C)<tr{gC). 

The other way around, a state a is separable, if and 
only if for aU C: fAB{C) > tr(d-C'). This is a kind of 
distance criterion. Our entanglement theorem [1] does no 
longer require the explicit form of any entanglement wit- 
ness. Entanglement is completely verified by Hermitian 
operators C & P, which are given by Eq. ([4]). The needed 
functions fAB{C) are readily obtained from Eq. ([5]). 

Let us now consider a bounded Hermitian operator A, 
which can always be expressed in terms of a positive op- 
erator C and a real number k e R, 

A^kI + C. (6) 

It is obvious, that all bounded Hermitian operators can 
be written in the form of A. This can be used to further 
simplify the theorem [T] 



3 



Theorem 2 A state g is entangled, if and only if there 
exists a Hermitian operators A: fAsi-^) < tT{gA). 

Proof. Note that tr(^) = 1. The function fAB is 
fAsiA) = sup{«; + (a, b\C\a, b) : {a\a) = {b\b) = 1} 

= K + fAB{C). 

From theorem]^ follows: 

fAB{C) < tr{ gC) 
n + fAB{C) <n + tY{gC) 
fAB{A) < trig A). 

From theorem [2] entanglement can be numerically 
tested. The set V and the set of bounded Hermitian 
operators have the same cardinality. Now the construc- 
tion of positive Hermitian operators, Eq. ([4]), becomes 
superfluous and the number of tests does not increase. 

Let us consider a simple implication of Theorem[2l The 
entanglement test for a bounded Hermitian operator A 
reads as 

sup{(a,6|i|a,6) : (a|a) = {b\b) = 1} < tr(^i). (7) 

The entanglement condition for the state g with the op- 
erator —A is 

snp{{a,b\[-A]\a,b)} < tr(^[-i]) (8) 
<^ M{{a,b\A\a,b)}>tT{gA). (9) 

Equation Q is a second entanglement condition for the 
operator A and it is equivalent to the original condi- 
tion ^ for the operator —A. 

Entanglement witnesses of the form W = A — 
inf{{a,b\A\a,b)}i had been considered before, see p^ . 
Here we gave the proof that any entanglement witness 
can be given in this form, which is a much stronger 
statement. This has been done for arbitrary dimensional 
Hilbert spaces. 



III. SEPARABILITY EIGENVALUE PROBLEM 

Let us consider the calculation of the function fAB- 
The Hermitian operator A has the following projections: 

Aa ^ tr a[A {\a) {a\ ® i b)] = {a\A\a), (10) 

ib ^tYB[A{iA'E)\b){b\)]^{b\A\b). (11) 

Now the extrema of {a,b\A\a,b) can be obtained. From 
Eq. (O, the extremum of the function 

gi\a), {a\, {b\) = g{a, b) = {a, b\A\a, b) (12) 

is calculated under the constraints hi{a) = {a\a) — 1 = 
and /i2(6) = {b\b) — 1 = 0. The functional derivatives are 



denoted by and This leads to two Lagrange 

multipliers Li,L2 and the conditions 



dg dhi dh2 



d{a\ 



-d{a\ 



■d{a\ 



(13) 



dg dhi 9/i2 



This can be written as: 



Ab\a) =Li\a), 
Aa\b) =L2|fo). 



(15) 
(16) 



Multiplying the first equation with (a| and the second 
with (&|, the Lagrange multipliers are obtained as Li = 
L2 = g. This leads to 

Definition 1 The equations 

Ab\a) g\a), 
Aa\b)^g\b), 

with the constraints {a\a) = {b\b) = 1, are called separa- 
bility eigenvalue equations. 

The separability eigenvalue problem can be solved by 
computers or, in simple cases, by hand. The closed set 
{|a)(8)|6) : (a|a) = {b\b) = 1} is bounded in Ha The 
smooth function g(a, b) defined on this set is bounded by 
3(0,^)1 < II ^11- Therefore a solution of the separability 
eigenvalue equation exists. According to the calculation 
(fT^ ~ we can formulate the following 

Lemma 2 The function fAB can be written as 

fAB{A) = sup{.g}, 

with g being eigenvalues of the separability eigenvalue 
equations. 

Proof. See the calculations ilS^) - jib]) and Eq. (S^. 

We also obtain inf{(a, 6| A|a, 6)} = infj^}, with g sepa- 
rability eigenvalue. In the following we want to consider 
some properties of the solution of the separability eigen- 
value equations. 

Proposition 1 Let go, |ao,t'o) be a solution of the sepa- 
rability eigenvalue problem of the bounded Hermitian op- 
erator A. 

1. The vector A|ao,6o) has a Schmidt decomposition, 
see l^], with the term golaoi&o); 

A|ao,6o) = 5o|ao,6o) + ^ ^kAO'k,bi), 
{\o.k)}k, {\i'i)}i orthonormal bases 

2. If gi ^ go, |ai,6o) (or |ao,6i)j is another solution 
then {ao\ai) = (or {bo\bi) = 0). 
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3. If gi 7^ go, |ai,6i) is another solution then |ao,6o) 
and |ai,6i) are linearly independent. 

Proof. 1. The general form of the vector is {tp) = 
A\ao,bo) = J2k,i^k,i\ak,bi). 

{ao, 6o|^|ao, ^o) = go = V'o.o 
{ak,bo\A\ao,bo) = (afe|AJao) 

= go{ak\ao) = = ipkM (k ^ 0) 
{ao,bi\A\ao,bo) = {bi\Aao\bo) 

= go{bi\bo) = ^ iJo,i (1^0) 



separability eigenvalue equations can be written for each 
component as 



X! '^k,i0.kbi ^ b*'^lJ^J = gai (Vi), (18) 

k,l J J 

'^ki^kk 9bj (Vj). (19) 

k,l J i 

Inserting Eqs. (fT5|) and into each other and using 
Eq. (fT2l) in the form 



Thus, = go\ao,bo) + J2k^o.i^o'^k,i\ak,bi). 

2. The first part of separability eigenvalue equa- 
tions read as 

AboWo) =5o|ao) and Abg\ai) = gi\ai). 

These are eigenvalue equations for the Hermi- 
tian operator Abg - acting on Ha - for differ- 
ent eigenvalues. Thus, (ao|ai) = 0. 

3. Assume linear dependence, a\aQ,bo) = 
|ai, bi). Since (ai|ai) — {bi\bi) — 1, we obtain 
\a\ — 1. This would imply (ao, bo|^|cto, ^o) — 
{ai,bi\A\ai,bi) , which is a contradiction to 
91 7^ 90- 

These properties of the separability eigenvalue equations 
can be easily seen for the solution of the example given 
in Sec. WJ\ 



IV. IMPLEMENTATION OF ENTANGLEMENT 
TESTS 

In the following we want to study two aspects of the re- 
sults obtained so far. In Sec. IIV Al an analytical solution 
of the separability eigenvalue problem will be derived for 
a special class of projection operators. In Sec. IIVBI a 
general entanglement test for arbitrary quantum states 
is under study. 



(20) 



k,l 



for g ^ (g = being a trivial case) we can separate 
Eq. 'lIlD from Eq. ^ 



(21) 
(22) 



With the interpretation ^ M we get 

MM^\a)^g\a), (23) 
MtM|6) = g\b). (24) 

The positive and compact operators MM^ and M can 
be given in spectral decomposition as 



MMt = ^|m,na,)(aJ, 

9 



(25) 
(26) 



Thus the non-trivial separability eigenvalues are gq,q 
Using M = Y,qrnq\aq){bq\ and |a'^) 



joq}, the state reads as 



(27) 



A. Tests by pure states 

For example, let us solve the separability eigenvalue 
equations for the special class of operators of the form 
A^ — \^p){'>p\. The normalized vector lip) can be ex- 
panded as 

In the same way the vectors \a) G Ha and ]b) G Hb can 
be written as \a) = X^i'^ik*) \b) ~ J2jbj\fj)- The 



where \a'q) and \bq) are orthonormal in each Hilbert space. 
This is the Schmidt decomposition of \^), cf. e.g. [l|]. By 
the above calculations we get the solutions 



|m,|^ p^q \a'q,bq) 
PJ^q \a'bq) ' 



fAsiA^,) = sup{|m,|^}. 



(29) 



For fAB{A^) — 1 the state — \a,b) is factorized, and 
A^ does not detect any entanglement. In all other cases, 
A^ is useful to identify entanglement. 
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Now we can write the special condition for entangle- 
ment, by use of in theorem [21 as 



'|/5|V) > sup{\mq\^}. 



(30) 



Let us consider the example of two harmonic oscillators 
in a mixture of a superposition of coherent states, \x~) — 
JV{\a,(3) — I — a, — /?}), with vacuum, |0,0}, 



emix-r,|x-)(x-l + (l-^)|0,0)(0,0| 



(31) 



where < 77 < 1 and A/" [2(1 - e-2(l"l'+l^'l'))]-i/2. 
For this state, higher-order moments are needed to verify 
NPT-entanglement even for r; = 1, see Based on the 
Bell state \^) = |$) = ^(|0, 1) + |1,0)), with /^s(i$) = 
1/2, the condition ([50)1 reads as 



V > 



sinhdaj 



(32) 



It identifies entanglement of Pmix for certain values of a, 
(3 and 77. For other choices of \^), even the simple con- 
dition pop may identify entanglement for more general 
parameters of the mixed state under study. 



with Apqrs = ^*spq- We obtain each element of this set 
{Ai\i=i,,,n for instance by varying each Ap^rs in the fol- 
lowing way. The absolute value |ApgVs| can be increased 

from to 1 with steps Ar". The argument arg(Apgrs) can 
be increased from to 27r with steps Ac/). Each oper- 
ator has the norm ||Ai||max < 1- We neglect operators 
with 1 1 Umax — and renormalize the other operators 
as Using the definition of the norm, we ob- 

tain that for each A with ||A||max — 1 exist one element 
Ai with 11^ - A^llmax < \/^r'^ + A(/)2 = e. In this way 
we may construct the test operators Ai for a desired pre- 
cision e. 

Now one can solve with an appropriate algorithm the 
separability eigenvalue equations for each Ai. This gives 
the values of Jab [Ai ) and according to equation ([33]) the 
optimal witnesses 



W,^ fAB{Ai)\-A, 



(38) 



Let us consider a grid of only 6 Hcrmitian operators. Fig- 
ure 1 indicates, to which extend these optimal witnesses 
identify entanglement. The gray area represents those 
entangled states, which are not identified. 



B. General test operators 

Let us now deal with the general form of theorem[2l We 
may explicitly construct all optimal entanglement wit- 
nesses as 

W^opt -/AB(i)i-A (33) 

More generally, any element of Was can be written as 
= Ai -I- Wopt (A > 0). Note that, if 7 is a positive 
real number and W an entanglement witness, then 
is also an entanglement witness, which is as fine as W . 

In the following we consider an implementation of our 
method for a finite dimensional Hilbert space T-La <S)'Hb- 
A numerical implementation could be a set of Hermitian 
operators {v4i}i=i...„ with the properties: 

lli.ll-l, (34) 
Vi: ||i|| = 1 3i ; <e. (35) 

This is a kind of a spherical grid, with a maximal dis- 
tance e > for any Hermitian operator A to at least one 
element Ai, for arbitrarily small e. 

For example let us consider the following construction. 
We use the most general form of a Hermitian operator A 
together with the norm || • Umax, 

^= ^pgrsl^p, fq){£r, fs\, (36) 

max — max l^^pg^sl, (37) 




FIG. 1: The chosen optimal witnesses Wi identify a manyfold 
of entangled states, except those in the gray areas. The set 
of separable states is approximated by the Wi. 



The test for entanglement is connected with an error 
depending on the chosen value of e. If e becomes smaller, 
the gray area in Fig. 1 is decreasing. For any given entan- 
gled state g there exists an e > 0, so that its entanglement 
can be identified by at least one of the chosen Hermitian 
operators Ai. In practice, the possible values of e can 
be related to the available experimental precision. The 
construction outlined above shows clearly that we need 
only a finite number of entanglement tests for the desired 
precission. 
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V. VERIFYING BOUND ENTANGLEMENT 

As we have mentioned above, for NPT states all wit- 
nesses can be given as C^'^ . Let us define the function 
gPT{a,b) = {a,b\C^'-^\a,b). The following proposition 
shows that the solution of the separability eigenvalue 
problem for C^"^ becomes superfluous. 

Proposition 2 If a general bounded Hermitian operator 
A has the solution g, \a,b) of the separability eigenvalue 
equations, then A^'^ has the solution g, \a, b*). It follows 

fAB{A^^)^fABiA). 

Proof. Since for the Hermitian operators \b){b\'^ = 
we find for all g(a, b) = (a, 6|A|a, b) 

g{a,b)=tY[A{\a){a\^\b){b\)] 

^tY[A^^ {\a){a\®\b*){b*\)] 
= gpT{a,b*). 

Thus, the optimization will lead to the solutions 
g, \a,b*). 

We find for the operators A^i, — \il'){tp\ presented in 
Sec. IIV Al the function fAsiA^^) = sup{\mq\'^}, with 
\mg\ being the Schmidt coefficients. Since for this exam- 
ple 5 = is a separability eigenvalue as well, see Eq. (p8)) . 
we get 

M{{a,b\A^\a,b)} = M{{a,b\Al^\a,b)} = 0. (39) 

To obtain the known characterization for NPT en- 
tangled states, gNPT, we use the entanglement con- 
dition as presented in Eq. which simplifies to 
tr(^NPT[|V')(V'r^)<0. 



Note that an arbitrary operator A^^ is a Hermitian 
operator as well. Since we can shift any operator A + nl 
without changing the entanglement witness, 

W-opt = fABiA)i - i = fABilii + i)i - (acI + i), (40) 

any entanglement test can be performed with an operator 
of the form C^'^. All entangled states which remain non- 
negative under FT are bound entangled ones, see p^ . 
The following characterization of PPT bound entangled 
states can be given: 

^BE is PPT bound entangled 

(1) VC'e-P: < tr(^BEC'^^) 

(2) 3C eV: inf{(a, b\C\a, b)} > tr(^BEC'^'^). (41) 

The first condition refers to PPT. The second condition 
identifies entanglement. The difference between C^"^ as 
a witness for NPT entanglement and for entanglement in 



general is equal to inf{(a, 6|C|a, 6)} - the minimal sepa- 
rability eigenvalue, see Fig. 2. 




FIG. 2: The gray area is the set of PPT bound entangled 
states. From the separable side, a typical witness is W = 
C^'^ - mi{{a,b\C^'^\a,b)}i. From the NPT side, a typical 
witness is IFnpt ~ (7^^. 



VI. SUMMARY AND CONCLUSIONS 

In the present paper we have proven the general form of 
entanglement witnesses. On this basis we have derived 
necessary and sufficient conditions for bipartite entan- 
glement. Optimal entanglement inequalities have been 
given in the most general form. They have been formu- 
lated with arbitrary Hermitian operators, which are easy 
to handle because of their well known structure and they 
are useful for applications in experiments. 

Separability eigenvalue equations have been formu- 
lated. They serve for the optimization of the entan- 
glement conditions for all chosen Hermitian operators. 
Some properties of the solution of these equations have 
been analyzed. The separability eigenvalue problem re- 
sembles the ordinary eigenvalue problem of Hermitian 
operators, with the additional restriction that the solu- 
tion is a factorizable vector. 

We have analytically solved the separability eigenvalue 
equations for a special class of projection operators. Us- 
ing these solutions, we could demonstrate entanglement 
of a mixed state given in terms of continuous variables. 
A general entanglement test of the proposed form can be 
implemented numerically, with its error being related to 
the available experimental precision. 

The method under study can also be used to identify 
any bound entangled state with a positive partial trans- 
position. This requires to test the given states for the 
negativity of its partial transposition. It turned out that 
the separability eigenvalues remain unchanged under par- 
tial transposition. Eventually, bound entanglement can 
be demonstrated by a combination of a general entan- 
glement test and a test for the negativity of the partial 
transposition of the state under study. 
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